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ABSTRACT 

Observations of gravitational microlensing in multiply imaged quasars currently pro- 
vide the only direct probe of quasar emission region structure on sub-microarcsecond 
scales. Analyses of microlensing variability are observationally expensive, requiring 
long-term monitoring of lensed systems. Here we demonstrate a technique for con- 
straining the size of the quasar continuum emission region as a function of wavelength 
using single-epoch multi-band imaging. We have obtained images of the lensed quasar 
MG 0414+0534 in five wavelength bands using the Magellan 6.5-metre Baade telescope 
at Las Campanas Observatory, Chile. These data, in combination with two existing 
epochs of Hubble Space Telescope data, are used to model the size of the continuum 
emission region a as a power-law in wavelength, <r a A". We place an upper limit on 

the Gaussian width of the r'-band emission region of 1.80 x 10 16 /if 1/2 ((M) /Me) 1 ' 2 cm, 
and constrain the power-law index to 0.77 ^ v ^ 2.67 (95 per cent confidence range). 
These results can be used to constrain models of quasar accretion discs. As a example, 
we find that the accretion disc in MG 0414+0534 is statistically consistent with a 
Shakura-Sunyaev thin disc model. 
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1 INTRODUCTION 

Models of quasar accretion a re currently poorly constrained 
by observations (Blaes 2007). Accretion discs are much too 
small for direct imaging, and so other observational meth- 
ods are required to probe their structure. Gravitational mi- 
crolensing of multiply imaged quasars is one such technique, 
uniquely able to probe accretion disc structure on scales 
~ 10 16 cm. 

The quasar central engine is thought to consist of 
an accretion disc surrounding a supermassive black hole. 
It is a general feature of accretion disc models (e.g. 
IShakura fe Sunvaevl Il973l : iHubenv fe Hubenvl Il997l ) that 
longer wavelength radiation is emitted from larger regions 
in the disc. Accretion discs are thus a perfect target for 
gravitational microlensing analyses, as the amplitude of fluc- 
tuations produced by gravitational microlensing depends 
strongly on the size of the emission region in the source. 



* This paper includes data gathered with the 6.5 meter Magellan 
Telescopes located at Las Campanas Observatory, Chile, 
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In particular, microlensing mod els predict larger fluctua- 



tions at shorter wavelengths (e.g. Rauch & Blandford 1991 
Jaroszvriski. Wambsganss fe Paczvriskill 19921 : Agol & Krolik 



1999; Yonchar a et al 



1999). 



Considerable effort has been invested in microlens- 



inasar emission 



regions (e.g. Wambsganss. Paezvriski & Schneider 1990b; 


Rauch & Blandford! 


199ll: IWitt. Mao & Schechterl 


1995; 


Wvithe et alJ 12000b 


: IWavth. O'Dowd & Websterl 


2005; 


Morgan et all 12006: 


Chartas et al. 2008h Recent efforts 



have focussed on observations in (typically) two wavelength 
bands, in an effort to better constrain models of quasar ac- 
cretion. The se analyses frequ ently make use of the method 
presented in iKochanekl l|2004l ). 

iPoolev et all (|2007m analysed ten quadruply imaged 
quasars (including MG 0414+0534, the quasar studied in 
this paper), and argued qualitatively that optical emis- 
sion regions must be larger than expected from basic 
thin-disc models by factors of ~ 3 — 30. This result 
was supported by iMorgan et alJ l|2007l h in which accre- 
tion disc sizes were determined for ten lensed quasars, and 
I Morgan et al.l (|200ct ). which focussed on X- r ay an d optical 
monitoring of PG 1115+080. lAnguita etaLI (120081 ) used g'- 
and r'-band monitoring of Q2237+0305 to obtain a spec- 
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tral slope consistent with the standard thin-disc models. 
IPoindexter. Morgan fc Kochanekl (|2008T ) also found an ac- 
cretion disc scale and slope consistent with standard thin 
disc models, using 13 years of photometric data of the two- 
image lens HE 1104-1805. 

The quantitative analyses discussed above rely upon 
monitoring of microlensing-induced brightness fluctuations 
over an extended time period. Constraints on the size of 
the emission region in the source are obtained by produc- 
ing a large number of simulated microlensing light curves, 
and searching for fits to the observational data. These con- 
straints are therefore dependent on the trans verse velocity, 
which can be estimated fro m the light curve (IWvithe et all 
l2000al ; lKochanekll2004r ). In lBate. Webster fc Wvithei l|2007l . 
hereafter BWW07) we demonstrated an alternative tech- 
nique for placing constraints on the size of quasar emission 
regions using single observations of lensed quasars displaying 
a flux ratio anomaly. 

Several gravitationally lensed quasars are observed to 
contain a pair of images straddling a critical curve. Basic 
lensing theory predicts that these images should have the 
same magnification. However in several cases, one image 
is obs erved to be significantly demagnified relative to the 
other JChang fc Refsdall 11973 ; iBlandford fc Naravanl [l98rj ; 
I Witt et all 1 19951 ). It has been argued that microlensing by 
a combination of smooth and clumpy matter in the lensing 
;alaxy is the most likely expla nation for this discrepancy 
Schechter fc Wambsganssll2002l ). 

By exploring microlensing simulations with varying 
source size and smooth matter content in the lens, we 
have previously obtained a 95 per cent upper limit on the 
size of the /-band emission region in the anomalous lensed 
quasar MG 0414+0534 of 2.62 x lQ 16 h^' 2 (M/M©) 1/2 cm 
(BWW07). As with all lengths measured from microlensing 
simulations, this size is normalised to the Einstein Radius, 
which is scaled by the square root of the microlens mass M. 

In this paper, we extend our analysis across multiple 
filters. Images of MG 0414+0534 in the r i', z', J and H 
filters were obtained on 2007 November 3, using the Mag- 
ellan 6.5-metre Baade telescope. We also make use of Hub- 
ble Space Telescope (HST) data in the NICMOS F110W 
and F205W filters, taken on 1994 August l4E and the 
WFPC 2 F675W and F814W filters taken on 1994 Novem- 
ber 8 (|Falco. Lehar fc Shapirdll997r i. These data allow us 
to obtain emission region size constraints in each filter or, 
equivalently, a power-law fit to the radius of the continuum 
emission region in the quasar as a function of emitted wave- 
length. 

The observations are described in Section 2. In Section 
3 we summarise our simulation method, with particular em- 
phasis on how it varies from that presented in BWW07. In 
Section 4 we present the resulting power-law fit to accretion 
disc radius as a function of emitted wavelength, and discuss 
our conclusions. 

Throughout this paper we use a cosmology with Ho = 
70kms" 1 Mpc~ 1 , Q, m = 0.3 and Ov = 0.7. 



Table 1. MG 0414+0534 magnitudes 



Filter 



source 


H 


J 


z' 


i> 


r' 


Ai 


17.87 


19.41 


21.92 


23.41 


25.40 




17.44 


18.85 


20.76 


21.97 


23.72 


B 


20.31 


21.51 


27.20 


28.71 


30.01 


C 


21.12 


22.84 


28.37 


29.51 


30.87 


L 


19.27 


20.59 


25.91 


27.72 


31.08 



Apparent magnitudes of the lensed images (Ax, A2, B, and C) 
and the lensing galaxy (L) in MG 0414+0534. Observations were 
taken in five filters with the IMACS and PANIC instruments on 
the Magellan 6.5-m Baade telescope, 2007 November 3. 



2 OBSERVATIONS 

We observed MG 0414+0534 using the Magellan 6.5-metre 
Baade telescope on 2007 November 3, in < 0'.'5 seeing. Ob- 
servations were made using the IMACS (in f/2 mode) r' , i', 
and z' filters and the PANIC J and H filters. Integration 
times were approximately five minutes per filter. Images are 
provided in Figure [1] 

We fit Gaussians to each quasar image. We began by 
assuming a round, Gaussian PSF, and then allow the ellip- 
ticity to vary as a free parameter. The //-band image best 
shows the separation of the lensed images, and this was ini- 
tially used to constrain the position of each component in 
each filter. In the final fit we freed up position as well, re- 
sulting in a 5-D fit (x, y, amplitude, FWHM, and ellipticity) 
for each component using a downhill gradient technique. 

The apparent magnitudes of the four lensed images and 
the lensing galaxy observed on 2007 November 3 are pro- 
vided in Table [T] Flux ratios in each filter between the two 
images of interest, F ts = A2/A1, are provided in Table [2] 
This table contains flux ratios from the observations pre- 
sented in this paper, as well as two sets of HST observa- 
tions, obtained fr om the CASTLES Survey web page 1 and 
iFalco et all l|l997f ). 

Table[2]shows that the flux ratio between images A2 and 
Ai at any single epoch grows larger as wavelength increases. 
We have assumed that this anomaly is solely due to mi- 
crolensing. An alternative explanation is differential extinc- 
tion due to dust in the lens. It has been argued previously 
that dust reddening is able to complet ely explain the anoma- 



lous flux ratio s in MG 0414+0534 dLawrence et al 



iMcLeod et"ail ll99Sl; lOstman. Goobar fc Mortselll boOcf ) 



1995; 



http: / /cfa- www. harvard.edu/castles/ 



However, iTonrv fc Kochanekl (|1999l ) showed that the lens 
galaxy is a passively evolving early-type galaxy, and thus 
unlikely to be dust dominated. Furthermore, several lens- 
ing systems displaying a flux ratio anomaly similar to MG 
0414+0534 are known to exist. Microlensing models predict, 
as is observed, that in such systems it is always the saddle 
point (A2) image that will be suppressed. Differential ex- 
tinction due to dust would not be expected to favour one 
image over another. 

Broadband observations of quasar continuum emis- 
sion are likely to be contaminated by broad emission lines 
(BELs), which may in turn affect a power-law fit to the ra- 
dius of the accretion disc. As we expect the continuum and 
broad emission line regions to be differently microlensed, 
this contamination is difficult to quantify exactly. Using 



The accretion disc in MG 0414+0534 3 




Figure 1. Magellan IMACS and PANIC imaging of MG 0414+0534, taken on 2007 November 3. The anomalous flux ratio F ob3 increases 
as we move blueward. 



Table 2. MG 0414+0534 flux ratios 



Filter 


Central A (A) 


fobs A2 


Date 


H 


16500 


0.67 + 0.05 


2007 November 3 


J 


12500 


0.60 ± 0.2 


2007 November 3 


z' 


9134 


0.34 + 0.1 


2007 November 3 


%' 


7625 


0.26 ± 0.1 


2007 November 3 


r> 


6231 


0.21 ± 0.1 


2007 November 3 


F110W 


11250 


0.64 + 0.04 


1997 August 14 


F205W 


20650 


0.83 + 0.03 


1997 August 14 


F675W 


6714 


0.40 + 0.01 


1994 November 8 


F814W 


7940 


0.47 + 0.01 


1994 November 8 



Central wavelengths and observed (anomalous) flux ratios F a i, s 
between images A2 and Ai in each of nine filters. The 2007 
November 3 observations were taken with the IMACS and PANIC 
instruments on the Magellan 6.5-m Baade telescope. The 1997 
August 14 observations were taken with the NICMOS instru- 
ment on HST (obtained from the CASTLES Survey web page 
http:/ /cfa-www. harvard.edu/castlcs/ ). The 1994 November 8 ob- 
servations w ere taken with the WFPC2 instrument on HST 
llFalco et al.lll9971) . 



the SP SS composite quasar spectrum of lVanden Berk et al.l 
(2001), we determined the BEL flux contribution to a 
(generic, unlensed) quasar at the same redshift as MG 
0414+0534 in our wavebands to be less than 5 per cent. 
The composite spectrum is shown in Figure [2] with the fil- 
ter transmission curves overlaid. The top panel shows the 



five Magellan IMACS and PANIC filters, the middle panel 
shows the two HST NICMOS filters, and the bottom panel 
shows the two HST WFPC2 filters used in our analysis. 

Optical and IR spectra are available for MG 
0414+0534 jHewitt et all 1 19921 ; lAngonin-Willaime et ail 
1994; L awrence et al.lll995l ). In these spectra, no BELs are 
detected within our wavebands. We therefore use the com- 
posite spectrum to illustrate the approximate size of BEL 
contamination. A detailed analysis of this source of error 
would require spectroscopy taken contemporaneously with 
each imaging epoch. 



3 SIMULATIONS 

We parameterise the frequency-dependent physical size of 
the accretion disc using Gaussian profiles of width a as a 
function of emitted wavelength A: 



a = ao 



(£)" 



(i) 



The observational data provide us with flux ratio as a func- 
tion of wavelength, F(X), at three separate epochs. Mi- 
crolensing simulations, detailed below, offer the link between 
source size and flux ratio, F(a). 

We conducted our microlensi ng simulations using an in- 
verse ray-shooting technique (e.g. Kavser. Refsdal fc Stabell 
ll986j ; IWambsganss. Paczvnski fc Katzlll990al ). The key pa- 
rameters for such simulations are the convergence Ktot and 
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Table 3. Lcnsing parameters 



10" 2x10" 
redshifted wavelength (A) 



3x10" 



Figure 2. Composite spectrum constructed from 2204 
quasar spectra selected from the Sloan Digital Sky Survey 
llVanden Berk et al.ll20o"lT) . Overlaid (dashed line) are the trans- 
mission curves for the IMACS r' , i' , z' filters and the PANIC J 
and H filters (top panel), the HST NICMOS Ff 10W and F205W 
filters (middle panel), and the HST WFPC2 F675W and F8I4W 
filters (bottom panel). The transmission curves have been rescaled 
by a constant factor for ease of comparison. 



the shear 7 of the lens at the i mage positi o ns. W e used the 
MG 0414+0534 lens model of I Witt et all l|l995h . provided 
in Table [3] The convergence of the lens can be split into two 
components - a compact stellar component k* and a con- 
tinuously distributed component k c . We allowed the smooth 
matter percentage s — n c / ' utot in the lens to vary between 
per cent and 99 per cent, in 10 per cent increments. 

The microlenses were drawn from a Salpeter mass func- 
tion dN/dM oc M" 2 - 35 with M max /M min = 50. Physi- 
cal sizes are therefore scaled by the average Einstein Ra- 
dius projected on to the source plane r)o, which is 3.75 x 
10 16 h TO l/2 ((A/)/M Q ) 1/2 cm for MG 0414+0534. Magnifica- 
tion maps were generated covering an area of 24r/o x 24r/o, 
with a resolution of 2048 x 2048 pixels. Twenty maps were 
generated for each image and for each smooth matter per- 
centage. 

We then randomly selected source positions in each of 
the Ai and A2 magnification maps. At each source position, 
we calculated the magnifications of 40 Gaussian sources of 
increasing characteristic width a. The characteristic width 
of the sources varied from 0.05?7o to 2.OO770 in increments 
of O.O5770. Dividing Ai magnifications by A2 magnifications 
allowed us to build up a library of 25500 curves of magnifi- 
cation ratio (equivalently, flux ratio) as a function of source 
width, Fi(a). 



Image 




7 


Mtot 


Ai 


0.472 


0.488 


24.2 


A 2 


0.485 


0.550 


-26.8 



Lensing paramet ers for images A\ and A2 in MG 0414+0534 
jWitt et al.lll995l) . 



We allowed the width of the inner-most (r'-band) emis- 
sion region ao to vary between 0.05?io and 2.00r/o, and the 
power-law index v to vary between and 15. These values 
were found to cover the interesting areas of parameter space 
in coarse simulations. For each combination of ao and v, we 
calculated emission region widths a„ at the central wave- 
length of each filter, A n (see Table [2}, using Equation [T] The 
value of Ao was taken to be 6231 A, the central wavelength 
of the r'-filter. 

By comparing a n (X n ) with each of our library of Fi(a) 
curves (interpolating linearly in a), we obtained 25500 val- 
ues of F™(A„) for each combination of ao and v at each 
epoch. These were then compared with F£ hB (\ n ), to obtain 
Xi values, which were converted to likelihoods using: 



Li(F° 



\ao,v,s) 



exp 



-Xi 
2 



(2) 



A final likelihood was obtained for each combination of ao 
and v at each epoch by summing over the individual likeli- 
hoods for each of the 25500 simulated flux ratio curves: 



L{F° bs \a ,v,s) = HF obs \a ,v,s). 



(3) 



We convert these likelihoods to an a-posteriori differen- 
tial probability distribution using Bayes' theorem: 



daodvds 



T / r^obs I \ dPp r i or 

oc L(F \ao,v,s)— 

dao 



dP dP 

prior ti - f p' 



dv 



ds 



(4) 



We assume constant Bayesian priors for inner width ao, 
power-law index v and smooth matter percentage s. 

Marginalising over smooth matter percentage, we ob- 
tain probability distributions for ao and v at each epoch: 



d 2 P 
daodv 



d A P 
daodvds 



-ds 



(•>) 



The individual probability distributions are provided in Fig- 
ure|3] The left panel shows the Magellan IMACS and PANIC 
distribution, the middle panel shows the HST NICMOS dis- 
tribution, and the right panel shows the HST WFPC2 dis- 
tribution. 

The probability distributions at each of the three epochs 
were then combined to obtain a final distribution for ao and 
v, which is shown in Figure U (top panel). 



4 RESULTS 

Figure U (top panel) shows contours through the prob- 
ability distribution for the Gaussian width of the inner- 
most (r'-band) emission region ao and power-law index 
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5 10 5 10 5 10 

v 

Figure 3. Probability distributions for width of the r'-band emission region ctq, in units of the average Einstein Radius T70 , and the 
power-law index v. Contours are 64, 26 and 14 per cent of peak probability. Distributions are shown for the Magellan IMACS and PANIC 
observations (left panel), the HST NICMOS observations (middle panel), and the HST WFPC2 observations (right panel). 



v. By marginalising over the power-law index v we ob- 
tain a 95 per cent upper limit on the width of the r'- 
band emission region in MG 0414+0534 of 0. 48770 (1.80 x 
10 16 /i7 1/2 ((M)/M©) 1/2 cm in physical units). Similarly, we 
can marginalise over 00 and obtain 95 per cent confi- 
dence limits on the power-law index of 0.77 ^ v ^ 2.67. 
Marginalised cumulative probability distributions are pro- 
vided in Figure |3] for <to (bottom left panel) and v (bottom 
right panel). The smooth matter percentage in this lens is 
unconstrained. 

The standard Shakura-Sunyaev thin-disc model has 
R oc A 4 / 3 (i gnoring the central temperature depression, 
which is located at scales smaller than smallest width probed 
in our analysis). This is consistent with our analysis at the 
68 per cent confidence level. 

Using a black hole mass for MG 0414+0534 of M bh = 
1.83 x 10 9 M q, obtained using the virial technique with 
the Hfi line I Peng et alj I2006T ) . we can calculate the pre- 
dicted half-light radius of the r'-band emission region from 
a Shakura-Sunyaev thin disc model. We can compare this 
theoretical value with our observational results as microlens- 
ing fluctuations are s ensitive only to the half-light radius of 
the e mission region l|Mortonson, Schechter fc Wambsganss! 
120051 ). We assume a 20 per cent accretion efficiency, and 
an accretion disc inner radius equal to the innermost stable 
orbit. We find that the Shakura-Sunyaev thin disc for MG 
0414+0534 in the r' band (1712A in the quasar rest frame) 
has a half-light radius R1/2 = 0.20r;o- This corresponds to a 
Gaussian source with characteristic width <to = 0.1 6770- This 
size can be combined with v = 4/3, and is marked as a cross 
in Figure |3] (top panel). 

There has been some disagreement between previous 
microlensing analyses regarding the size of o p tical emis- 
sion regions i n len sed quasars. iPoolev et all (|2007l ) and 
iMorgan et al.l (|2008l ) both found that quasar optical emis- 
sion regions are larger than pre dicted by Eddington lim- 
ited thin disc models. Conversely, lAnguita et al.l l|2008l ) and 



iPoindexter et all (|2008l ) found disc sizes consistent with thin 
disc models. 

IPoindexter et all (|2008l ) noted that the solution to this 
discrepancy might be a slightly shallower temperature pro- 
file than predicted for a Shakura-Sunyaev thin disc. Our 
results are consistent with this suggestion. We note, how- 
ever, that the onl y previous analysis of MG 0414+0534 
ijPoolev et al.| [2007) found that it was one of only two lensed 
quasars in a sample of ten for which a Shakura-Sunyaev disc 
was just large enough to explain both the observed X-ray 
and optical flux ratios . 

IPoindexter et all l|2008l ) also fit a power-law in wave- 
length to accretion disc radius for the lensed quasar HE 
1104-1805. They obtained a power-law index of v = 
l-64lQ'4g (68 per cent confidence). At the 68 per cent con- 
fidence level, the power-law index obtained here is v = 
1.481q'43. These results are indistinguishable at the la level. 
This is particularly interesting as the methods used in the 
two analyses differ significantly. IPoindexter et all (|2008l ) 
compared 13 years of monitoring data in 11 wavebands to 
simulated microlensing light curves. We take advantage of 
the fact that two images straddling a critical curve each 
behave differently when microlensed, allowing us to obtain 
similar results with only 3 epochs of imaging in 10 overlap- 
ping wavebands. 



5 CONCLUSIONS 

In conclusion, we have demonstrated a technique for us- 
ing near-contemporaneous multi-band imaging of anomalous 
lensed quasars to constrain the radius of the quasar contin- 
uum emission region as a function of emitted wavelength. 
The application of this technique to MG 0414+0534 yielded 
an upper limit on the width of the r'-band emission region 
of 1.80 x lO 16 /i7 O 1/2 ((M)/M ) 1/2 cm, and a power-law index 
range of 0.77 ^ v ^ 2.67, both with 95 per cent confidence. 
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Figure 4. Top: probability distribution for width of the inner- 
most (r'-band) emission region <rrj, in units of the Einstein Ra- 
dius rjo, and the power-law index v. Contours are 64, 26 and 14 
per cent of the peak probability. The cross marks the locus of a 
Shakura-Sunyaev disc with 20 per cent accretion efficiency and 
M bh = 1.82 X 10 9 M Q jPene et al.ll2006h . Bottom left: cumula- 
tive probability distribution for the width of the r'-band emis- 
sion region <ro, in units of the Einstein Radius rjo- Bottom right: 
cumulative probability distribution for the power-law index v. 



These results are consistent with a Shakura-Sunyaev accre- 
tion disc model. 

The greatest strength of our method is its observational 
efficiency. Previous analyses have relied upon multi-epoch 
monitoring, with the tightest size constraints typically ob- 
tained during the fortuitous observation of high magnifica- 
tion events. The method demonstrated here requires only 
as much time as is necessary to obtain sufficient signal to 
noise of the lensed images in each filter. However, it is only 
useful if the lensed quasar is displaying an anomalous flux 
ratio. We have assumed this anomalous flux ratio is pro- 
duced solely by microlensing, although we cannot rule out 
differential extinction due to dust in the lensing galaxy as 
an alternative explanation. 

We have obtained observations of additional systems 
with anomalous flux ratios, the analysis of which are forth- 
coming. 
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